Quantum Vacuum Friction in Highly Magnetized Neutron Stars 
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In this letter we calculate the energy loss of highly magnetized neutron star due to friction 
with quantum vacuum, namely Quantum Vacuum Friction (QVF). Taking into account one-loop 
corrections in the effective Heisenberg-Euler Lagrangian of the light-light interaction, we derive an 
analytic expression for QVF allowing us to consider magnetic field at the surface of the star as high 
as lO" T. In the case of magnetars with high magnetic field above the QED critical field, we show 
that the energy loss by QVF dominates the energy loss process. This has important consequences, in 
particular on the inferred value of the magnetic field. This also indicates the need for independent 
measurements of magnetic field, energy loss rate, and of the braking index to fully characterize 
magnetars. 
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INTRODUCTION 

Quantum Vacuum Friction (QVF) is a very fundamen- 
tal phenomenon [l[ related to the fact that quantum vac- 
uum can be regarded as a standard medium whith its 
own energy density and electromagnetic properties. Well 
known effects due to static vacuum properties are the 
Casimir effect J2| and the vacuum magneto-electric opti- 
cal properties [3|] . Recently, dynamical effects due to vac- 
uum viscosity have been studied [H , in particular photon 
radiation stimulated by moving mirrors [3] (i-C dynami- 
cal Casimir effect) and friction betweenperfectly smooth 
surfaces moving relative to each other 0, @, Q • 

Considering a classical magnetic dipole moment rotat- 
ing in a standard medium with magneto-electric optical 
properties, relaxation effects (at molecular level for ex- 
ample) produce retardation effects between the induced 
magnetization of the medium and the rotation of the 
magnetic dipole moment. In this case, the temporal delay 
between the medium's response and the inductor gives 
rise to a classical dissipative frictional force. 

In this letter, we investigate the friction effect between 
a non stationary magnetic dipole moment and its in- 
duced quantum vacuum magnetic dipole moment. In 
this case the temporal delay between the vacuum's re- 
sponse and the inductor is due to the finite velocity of 
light. We calculate the friction energy loss rate in the 
case of highly magnetized neutron stars. Pulsars are the 
more appropriate systems to look for such an effect. Pul- 
sars are fast rotating neutrons stars, with a very high 
magnetic dipole moment tilted with respect to their ro- 



tational axis 8] . Typically, their mass is of the order of 
that of the Sun, and their radius of the order of 10 km. 
The magnetic field of neutron stars is typically of the 
order of 10* — 10^ T, while in the case of magnetars |^ 
astrophysical observations seem to indicate that fields as 
high as 10"'^^ T exist on their surface 13, U, 12|. Spin- 
ning periods of neutron stars range from tens of millisec- 
onds for young stars to seconds, increasing with star age. 
Because of their fast rotation retardation effects are sig- 
nificant. A rotating neutron star can thus be considered 
as a classical magnetic dipole moment weakly coupled to 
a bath, namely the quantum vacuum described by QED 
as electron-positron pairs due to quantum fluctuations. 
In a microscopic description, the energy is dissipated via 
the polarization process and through the annihilation of 
the polarized pairs which relax this excess of energy to 
the bath. This excess of energy corresponds to the fric- 
tionnal energy we calculate using a classical interaction 
between the star's dipole moment and the bath. We show 
that, for neutron star with a very high magnetic field, 
that exceeds the Quantum ElectroDynamics (QED) crit- 
ical field {Be — 4.4 10^ T), the energy loss is essentially 
due to QVF, while the energy loss due to the classical 
process as star rotating dipole radiation becomes neg- 
ligible. This applies in particular to magnetars which 
associate a strong magnetic field and a few seconds spin- 
ning period. To infer the magnetic field of such stars one 
cannot use the classical dipole energy loss formula. Our 
results suggest that to characterize magnetars one needs 
to measure the energy loss rate and the magnetic field 
independently. Taking into account QVF, we calculate 
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the braking index [l3[ and show that its measurement 
can provide a non-model dependent determination of the 
magnetic field on the sm'face of a neutron star. 



ENERGY LOSS RATE 

To study the quantum vacuum magnetization in the 
presence of a magnetic field, we start with the effective 
Heisenberg-Euler Lagrangian Che of the light-light in- 
teraction Due to relativistic invariance, Hjje is a 
function of the Lorentz invariants, F and G, given by [l5| : 
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where ep is the vacuum permittivity and fio the vacuum 
permeability. When one-loop corrections are included, 
Che can be written as Che = Cq + Ci, where Cq ^ F/2 
is the usual Maxwell's term. Neglecting the Electric field, 
F = — B^//io and G = 0. For Ci we use the analytic 
expression derived by Heyl and Hernquist 1^, 17 1: 
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4.4 10^ T is the critical field and 
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with In A = 0.248 754 477. Denoting by Mqv the induced 
magnetization, (i.e the quantum vacuum magnetic dipole 
moment per volume element), we have [l8| : 
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Using the analytic form ^ of the Lagrangian, we can 
calculate Mqv to the first order of a = /he ~ 1/137: 
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FIG. 1: Rotating pulsar in vacuum. 6 is the inclinaison of the 
star's magnetic dipole moment with respect to its rotation 
axis, dr denotes a vacuum element of volume. 



Let us then now consider a neutron star rotating in vac- 
uum (see fig.[T|). We denote by m its magnetic dipole mo- 
ment, R its radius and Bq the magnitude of the magnetic 
field at its surface {Bq ~ ^pm/A-KR^ where m — |lm||). 
We define a fixed frame {x,y,z) with z-axis parallel to 
the rotation axis of this star. {0, ip) stand for the spher- 
ical polar angles of m in this fixed frame. If p denotes 
the spinning period of the neutron star, (p — cut, with 
Lu — 2tt/p and m.Uz mcos6. At time t the magnetic 
moment of the star produces a magnetic field B(r,t). Let 
us denote by r the position vector of a vacuum element 
of volume dr ~ sin f3drdf3d'y where (r,/3, 7) stand for 
the spherical coordinates of r in this fixed frame.. Since 
Lur/c <^ 1 inside the region where quantum vacuum mag- 
netization is important, the leading term within the dipo- 
lar magnetic field approximation is |19| : 
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In this expression, retardation effects have been taken 
into account with the argument t — r/c in m. Accord- 
ing to eq. ([5]), the induced quantum vacuum magnetic 
moment at r is given by : 
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At time t + r/c the magnetic field dBq^, produced by 
dmqy{r,t) at the center of the star is : 
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This field interacts with the magnetic dipole moment of 
the star. At this stage, quantum vacuum can be regarded 
as a standard medium. Therefore, energy loss rate due 
to this friction is given by the classical formula [3] ■ 

dEqy = - {m{t + r/c) x dBq„(0,< + r/c)) uj.u^ (12) 

The total QVF energy loss rate is obtained by integrating 
eq. over the space external to the star, and averaging 
it over time : 



r=R J 13=0 Jj=0 
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Eqy can thus be calculated numericaly. Since retarda- 
tion effects grow with r while the star's magnetic field 
decreases significantly, the space region that mainly con- 
tributes to QVF is located around few star's radii and in- 
side the light cylinder. Moreover, in the case of a dipolar 
magnetic field, using eq. one can show that for Bo < 
IQio T, iB^iv,t)/B^_)fqy{B^{v,t))^^ {AB^(v,t)/AhBl) 
and we obtain the analytic expression 
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FIG. 2: Quantum Vacuum Friction energy rate Eq^ (full line) 
and classical matter radiative energy loss rate Er (dotted line) 
as a function of Bo sinS for the Crab pulsar B0531-I-21. 

In fig. [2l we have plotted the function Eqy obtained 
numericaly from eq. (|13p versus Bq sin Q (full line) , for 
the Crab pulsar with a spinning period p = 33,11tos. 
We have also plotted the star rotation radiation rate Er 
versus Bq sin Q given by classical dipole model bdl (since 
in this case, QED corrections are not relevant [17|). 
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Fig. [2] clearly shows that for fields exceeding the QED 
critical field the star energy loss is essentially due to QVF. 

At this stage, it is important to stress that spinning 
period dependence is not the same for the classical case 
and for the QVF. The ratio of classical to QVF losses 
decreases like 1/p^ for large -p. QVF becomes more im- 
portant for slowly rotating neutron stars. 

This result has important consequences. Actually, the 
value of the magnetic field on the surface of a neutron 
star is inferred by the energy loss of the star, derived 
by the measured value of the spin-down rate, assuming 
that the energy loss of the star is given by classical dipole 
model. In a case for which QVF should play an impor- 
tant role, the inferred value of the magnetic field should 
then be modified. Moreover, we see that the method 
used to get the value of the magnetic field, gives Bq sin0 
rather than Bq. One should thus assume the value of 
B to get Bq (typically 9 is taken equal to 7r/2). In the 
following, taking advantage of the fact that this angular 
dependence is the same for Eqy and we show that 
the inferred value of the magnetic field at the surface of 
a pulsar can be derived by the measurement of the star's 
braking index. 



BRAKING INDEX 

The braking index n, is a fundamental parameter of 
pulsar electrodynamics describing the rate at which a 
magnetized neutron star loses rotational energy. This 
dimensionless quantity is given by 
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where v is the spinning frequency of the pulsar and z> (re- 
spectively i>), denotes the first (respectively the second) 
derivative of v with respect to time. From this definition 
one can see that the braking index can be determined 
from pulsar timing measurements without any assump- 
tion concerning the star's structure and so, provides a 
crucial information for our understanding of the physics 
underlying pulsar spin-down. From a theoretical point of 
view, the pulsar's rotational energy loss rate is given by 



(17) 



where / is the moment of inertia of the rotating star. 
Assuming a pure classical dipole energy loss mechanism, 
we obtain from cq. (|15p : 
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where Kr = Er/^n'^Iv'^- Taking the logarithmic deriva- 
tive of this equation, we find n = 3. So far, braking 
indices of only six pulsars have been precisely measured, 
all of which are remarkably smaller than the value n = 3 
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expected for pure classical dipole radiation model (see ta- 
ble U second column) . This result clearly shows that an 
additionnal energy loss mechanism should be included 
in the model to fit the data. Explanations of this dis- 
crepancy have already been suggested as for example the 
effect of the relativistic pulsar wind [2l| or the fact that 
the pulsar's magnetism can not be modeled with a sim- 
ple dipolar field [i^]- In this letter, we do not discuss 
the validity of such assumptions but we propose instead 
QVF as an additionnal energy loss mechanism which can 
provides a new coherente explanation for all the values 
of braking index measured so far. Taking into account 
QVF, eq. m]) becomes 
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where Kgy = Eq^/An'^Iv^. One can then derive the cor- 
responding braking index 
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which can be written as 
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From this last expression together with eq. (fT^ and 
eq. (fT5)) we see that the braking does not depend on 9. 
More precisely we get 



(22) 



In the case of magnetars [Bq 3> Be) we expect a value 
of n ~ 1. In fig. [2] we show how the inferred value of the 
magnetic field B'^^ can be obtained by matching the an- 
alytical expression given by eq. (|22p with the measured 
value of a pulsar's braking index. B^^^ is thus coher- 
ent with n and does not require any assumption on the 
pulsar's geometry. 

We have applied this method with the six braking in- 
dices wich have been measured with certainty. Results 
are presented table H] (third column) together with val- 
ues of i3g=^ the magnetic field obtained with a pure 
classical dipolar model assuming 6 = 7r/2 and n = 3 
(last column). Depending of the case, the inferred value 
Bq"'^ can be higher or smaller than Bdipoie- Neverthe- 
less, it is interesting to notice that all the values found 
for the inferred magnetic field are smaller than the value 
of the QED critical field Be, that is not the case with 
the classical dipolar model. This last result can provide 
a very important test to confirm QVF predictions. Be- 
cause of the gamma-ray emitting electron-positron cas- 
cades which occur for E above the threshold value of 
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FIG. 3: Derivation of the inferred value of the magnetic field 
at the surface of the Crab pulsar B0531-I-21 using the analyt- 
ical expression given by eq. (I22p . 



i? ~ 3 X lQt''^erg.s~^ 2^, lO], magnetars, i.e pulsars with 
a magnetic field above the QED critical field (within the 
classical dipolar model) should not be gamma emitters. 
Detecting a gamma-ray emission from such a system, for 
example by the NASA GLAST mission due to launch in 
2008, would provide an evidence in favor of QVF which 
leads to smaller values of the magnetic field for magnetars 
compatible with gamma-ray emission. 



Name n 


u 

is-') 


(10'^ G) (10" G) 


J1846-0258 2.65(1) 


3.07 


1.0 


49 


B0531-K21 2.51(1) 


30.2 


12.4 


3.8 


B1509-58 2.839(3) 


6.63 


1.4 


15 


J1119-6127 2.91(5) 


2.45 


0.4 


42 


B0540-69 2.140(9) 


19.8 


12.4 


5.1 


B0833-45 1.4(2) 


11.2 


16.1 


3.4 



TABLE I: Spin and inferred magnetic field f or p ulsars with 
measured braking index n taken from [2^. [2^, [251 [2^ |27|. [2^1 . 



CONCLUSION 

We have studied the energy loss due to friction with 
quantum vacuum in the presence of high magnetized neu- 
tron star, namely Quantum Vacuum Friction (QVF). 
In neutron stars with magnetic fields of the order of 
10^ — 10^ T, and spinning period of less than a second, 
this effect is negligible compare to classical energy loss 
dipolar process. In the case of magnetars with high mag- 
netic field, and spinning period of few seconds, we have 
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shown that the energy loss by QVF dominates the en- 
ergy loss process. This has important consequences, in 
particular on the inferred value of the magnetic field. It 
also indicates the need for independent measurements of 
magnetic field, energy loss rate, and of the braking in- 
dex to fully characterize magnetars. QVF should play as 
well an important role in the early stages of neutron star 
formation and it should Jdo taken into account in models 
of magnetar (see e.g. [31]), also for what concerns its 
spindown history. On the other hand, macroscopic viola- 
tion of the linearity of Maxwell's equations predicted by 
QED 3] has not yet been proved. Experimental search 
of vacuum birefringence is in progress (see (3^ and ref. 
within). Evidence for QVF in the neighbourhood of a 
neutron star would thus be the first confirmation of such 
fundamental QED prediction, while negative observation 
would be the first indication of limitations on the QED 
description of vacuum. 

The authors thank P. Jones, S. Sarkar and D. A. Smith 
for helpful discussions. 
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